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! Abstract 
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■ The abundance and size distribution of quark nuggets (QN), formed a few 

microseconds after the big bang due to first order QCD phase transition in 
| the early universe, has been estimated. It appears that stable QNs could be 

ON ■ a viable candidate for cosmological dark matter. The evolution of baryon 

inhomogeneity due to evaporated (unstable) QNs are also examined. 
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1 Introduction 

As per the standard model, the universe, after a few micro seconds of the big 
bang, underwent a phase transition from quarks and gluons to hadrons. There are 
well organised efforts to mimic such a phase transition in the laboratory through 
heavy ion collisions at ultra- relativistic energies [IJ. Although there are similarities 
between the two scenarios, the relatively long time scales in the early universe phase 
transition may indeed be more conducive to a reliable thermodynamic description. 
In the absence of any consensus on the order of QCD phase transition for two light 
and one medium-heavy quarks form lattice calculations || in the present work we 
assume an underlying picture of a first order phase transition |3j from quarks and 
gluons to hadrons. The central question we would like to address is what plausible 
1 Present Address: Physics Department, Kyoto University, Kyoto 606-8502, Japan 
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remnants may have survived since that primordial epoch. Crawford and Schramm M 
and Van Hove || argued that fluctuations in the horizon scale, triggered by the 
phase transition, may lead to the formation of primordial black holes, which could 
be as large as M , the solar mass. Schramm || has recently suggested that these 
black holes could even be the candidates for the Massive Compact Halo Objects 
(MACHO 's) |7|, ||, which had of late been discovered in the halo of the Milky way, 
in the direction of the Large Magellanic Cloud (LMC) using the gravitational lensing 
techniques. On the other hand, a first order QCD phase transition scenario involving 
bubble nucleation at a critical temperature T c ~ 100 — 200 MeV could lead to the 
formation of quark nuggets (QN), made of u, d and s quarks at a density somewhat 
larger than normal nuclear matter density. If these primordial QN's indeed survive 
till the present epoch, they could be a possible candidate for the baryonic component 
of the dark matter ||. Such a possibility would be aesthetically rather pleasing, as it 
would not require invoking any exotic physics nor would the success of the primordial 
(Big Bang) nucleosynthesis scenario be materially affected || [10], [EL], |l2]j . 

The central question in this context is, thus, whether the primordial QN's can 
be stable on a cosmological time scale. The first study on this issue was addressed 
by Alcock and Farhi [13] who argued that neutrons could be liberated and emit- 
ted from a QN at temperature T > I N , where 1^ (~ 20-80 MeV) is the neutron 
"binding energy" i.e. the difference between energy per baryon in strange mat- 
ter at zero temperature and the mass of the nucleon. They calculated the rate of 
baryon evaporation from QNs by using detailed balance between the two processes 
of neutron absorption and emission in a system consisting of QNs (of all possible 
baryon numbers) and neutrons. They also assumed a geometric cross section for 
the neutron absorption by the QNs and complete transparency of the nugget sur- 
face to neutrons. This gave the result that QNs even with the largest allowed initial 
baryon number (i.e. the total baryon number Ns^hor ~ 10 49 (100MeV/T) 3 contained 
in a horizon-size volume of the universe, at the time of formation of the nugget, T 
being the temperature of the universe at that time) were unstable with respect to 
baryon evaporation. These results apparently eliminated the possibility of any QN 
surviving till the present epoch. 



Madsen et al |14[ then pointed out that, since evaporation was a surface process, 



2 



emission of neutrons (proton emission was suppressed due to the Coulomb barrier) 
made the surface layer deficient in u and d quarks, although relatively enriched in 
s quarks. The rate of conversion of s quarks back to u and d quarks as well as 
convection of u and d quarks from the core of the nugget to the surface were both 
too slow to establish flavor chemical equilibrium between u, d and s quarks on the 
surface layer. As a result of this deficiency of u and d quarks in the surface layer, 
further nucleon evaporation was suppressed. Madsen et al found that QNs with 
initial baryon number Nb > 10 46 could well be stable against baryon evaporation. 
QNs at temperatures above a few MeV could also be subject to the process 



of "boiling" |]To[ , i.e., spontaneous nucleation of hadronic bubbles in the bulk of 
the nugget and consequent conversion of the nugget into nucleons. Madsen and 
Olesen [16| however, showed that although boiling is thermodynamically allowed, 
the time scale for bubble nucleation inside QNs is too long compared to the time 
scale of surface evaporation for reasonable values of the parameters used. 

All of the above studies used thermodynamic and binding energy arguments to 
calculate the baryon evaporation rate; the microscopic dynamics determining the 
probability of baryon formation and emission had been neglected. Clearly, for a 
realistic description of the process a dynamical model of baryon emission from QNs 
was needed. 



Earlier, Bhattacharjee et al [I7| used the chromoelectric flux tube model (in- 
spired by QCD) to demonstrate that the QN's would survive against baryon evap- 
oration, if the baryon number of the quark matter inside the nugget was larger 
than 10 42 . For reasons explained in ref. |T7|], this estimate is rather conservative. 
Sumiyoshi and Kajino |TJ| estimated within a similar approach that a QN with an 
initial baryon number ~ 10 39 would survive against baryon evaporation. 

It may also be mentioned in this context that Madsen and Riisager |12| had cal- 
culated the minimum size of the QNs from the constraint imposed by the primordial 
He abundance. The authors showed that the minimum radius of the QNs should 
be more than ~ 10~ 6 cm, otherwise the absorption of neutrons by the QNs will 
upset the neutron to proton ratio in the universe, resulting in a helium abundance 
in complete disagreement with the experimental results. 

In spite of these efforts, not much attention has been paid towards the issues 
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of formation and size-distribution of the surviving quark nuggets in the universe. 
The size-distribution and abundance of the QNs is very important in the context 
of their candidature as dark matter. The calculation of a lower cut-off in size, if 
any, would tell us what can be the minimum size and the baryon number content 
of a QN that we should look for. On the other hand, the distribution function also 
indicates the most probable size of the QNs. Such studies in the framework of the 
GUT (Grand Unified Theory) phase transition and the associated "Trapped False 
Vacuum Domains" (TFVD) have been done earlier . The basic contention of this 



paper is to carry out these studies for the QCD phase transition in the early universe. 
We also study the evolution of baryon inhomogeneity created by evaporating QNs, 
with baryon number lower than the survivability window, as mentioned above, due 
to the conduction of heat by neutrinos, and consequently the dissipation of the 
baryon inhomogeneities. 

We organise the paper as follows. In the next section we evaluate the size- 
distribution of the QNs. The results for different nucleation rates have been dis- 
cussed. Section 3 contains the evolution of baryon inhomogeneity which originated 
due to the unstable QNs. In section 4 we conclude. 

2 The size-distribution of quark nuggets 

The evolution of the universe during the QCD phase transition is governed by Ein- 
stein's equations, 

(R\ 2 _ 8tt P d(pR 3 ) dR* 

\RJ " 3mJ' dt dt 1 ' 

where p is the energy density, P the pressure and m p i the Planck mass. In the above 
equation, R is the cosmological scale factor in the Robert son- Walker space time and 
is defined by the relation 

ds 2 = -dt 2 + R 2 dx 2 = R 2 (-d£ 2 + dx 2 ) 

dx 2 = dX 2 + X 2 {sin 2 6d<j) 2 + d6 2 ) (2) 

where X is the co-ordinate radius i.e. the radius in the unit of the cosmological 
scale factor R(t). 
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It is well known that in a first order phase transition, the quark and the hadron 
phases co-exist in a mixed phase at the critical temperature of transition. Around 
the critical temperature, the universe consists of leptons, photons and the massless 
quarks, anti-quarks and gluons, described in our case by the MIT bag model with 
an effective degeneracy 51.25). The hadronic phase contains relativistic ir- 

mesons, photons and leptons with a small baryon content (ps/p 7 ~ 1CT 10 ) and 
is described by an equation of state corresponding to massless particles with an 
effective degeneracy = 17.25. The energy densities and the pressures of hadronic 
and quark matter are given by 
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where B is the bag constant. 

The evolution of the scale factor in the mixed phase is given by (see also f2~0| 
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In the process we also get the volume fraction of the quark matter f(t) in the mixed 
phase as 
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where r = p q / p^, t c = ^?>m 2 pl /d>ixB is the characteristic time scale for the QCD phase 
transition in the early universe and U is the time when phase transition starts. (The 
definition of r and hence the expressions of R(t) and f(t) here are somewhat different 
from that in ref. [^0|1 ). The characteristic time scale depends on the bag constant 
and hence on the critical temperature of the quark-hadron phase transition (T c ). In 
fact t c = 144/is for T c = lOOMeV and t c = 64/is for T c = 150MeV. 

In the coexisting phase, the temperature of the universe remains constant at T c , 
the cooling due to expansion being compensated by the liberation of the latent heat. 
In the usual picture of bubble nucleation in a first order phase transition scenario 
hadronic matter starts appearing as individual bubbles. With the progress of time, 
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more and more hadronic bubbles form, coalesce and eventually percolate to form an 
infinite network of hadronic matter which traps the quark matter phase into finite 
domains. The time when the percolation takes place is usually referred to as the 
percolation time t p , determined by a critical volume fraction f c , (f c = f(t p )) of the 
quark phase. 

In an ideal first order phase transition, the fraction of the high temperature 
phase decreases from the critical value f c , as these domains shrink. For the QCD 
phase transition, however, these domains could become QN's and as such, we may 
assume that the lifetime of the mixed phase tf ~ t p . 

As mentioned above, just after percolation one can have pockets of quark matter 
trapped as bubbles in the ambient hadronic matter. The probability that a spherical 
region of co-ordinate radius X at time t p with nucleation rate I(t) lies completely 



within the quark matter domain is given by |19 
P(X,t p 



cxp 



47T 



dtI{t)R\t) (X + X(t p ,t)Y 



(6) 



where X(t p ; t) is the coordinate radius of a bubble, at time t p , which was nucleated 
at time t. 

For convenience, let us now define a new set of variables z = XR(ti) / vt c , x = t/t c 
and r(x) = R(x) / R(xi); where v is the radial growth velocity of the nucleating 
bubbles. Then 



P(z, x p ) = exp 



— —vHt 



dxl{x) {zr(x) + y(x p , x))'' 



where 

rx 

y(x,xf) = / r(xf) /r(xff)dxff 
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So the fraction of quark matter present at time t p is 

r 47r 

f c = P(0,x p ) = exp 



„vt c dxl(x)y (x p ,x) 
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(7) 
(8) 

(9) 



Let us now look at the size-distribution of the Trapped Quark Matter Domain 
(TQMD). In order to do so we will follow the procedure of ref. fll9| . The difference 
of our work from that of Kodama et al. is that we have considered exactly spherical 
nuggets whereas they have included a deformation factor. It should however be 
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noted that the deformation factor, as found by Kodama et al, is small. Moreover, 
due to the presence of non-zero surface tension in case of QCD phase transition the 
bubbles are likely to be spherical. Even more importantly, we focus our attention to 
the percolation time t p when the hadronic matter forms the ambient background. 
All these considerations allow us to consider the false vacuum domains (the quark 



phase) as being spherical in shape. Following Ref. [19] let us assume that F(X; t)dX 
is the number of TQMDs per unit volume within the size {X, X + dX} at time t. 
Then P(X,t) can be thought to be the probability that a QN of coordinate radius 
X at a fixed position is contained in a TQMD. Now a TQMD of size 77 can contain 
such a sphere of size X only when the center of TQMD lies within the coordinate 
radius 7] — X from the center of the sphere. If a is the minimum size of a TQMD 
i.e. F(X,t) vanishes for X < a then one can write 

P(X;t)= —( V -X)*F(rr,t)dV (io) 

Ja+X 6 

The distribution function vanishes for X < a. One can now solve the above 
equation using Laplace transformation to obtain F(X) (see appendix). 
The result, in terms of z, is 
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The solution of the equation F(a) = gives the minimum size of the quark- 
nugget a. Now, the number of nuggets per unit volume is given by 



POO 

R- S (t p ) / F(X)dX 

J a 



= r-% p ) r^f^ (12) 

Jo V Z c 

The volume of each quark nugget is given by -^7i(zvt c ) 3 . Since the visible baryon 
constitutes only ten per cent of the closure density (Qb = 0.1 from standard big 
bang nucleosynthesis), a total of 10 50 baryons will close the universe baryonically 
at T = 100 MeV. We emphasize at this point that these QNs would not disturb 
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the standard primordial nucleosynthesis results to any considerable extent, as they 
would not participate in usual nuclear reactions. Therefore, if we assume that the 
total baryon content of the dark matter is carried by the quark nuggets, then, 

N B = 10 50 (100/T(Me\/)) 3 = V H ™™ p B ^ f{z)zHz (13) 

where Vh is the horizon volume and p B is the baryon density inside each nugget. 
We now solve the above equations self-consistently to obtain a, t p and f c . These 
values are then used to study the size-distribution of the quark nuggets. 

To calculate the size and distribution of QNs we need to know the rate of 
nucleation during the phase transition process. Many authors [23, 24" | have proposed 



various nucleation rates for the QCD phase transition. In the absence of a consensus 
as to which rate is preferable to others, we look at several of them. We begin with 
the rate proposed by Cottingham et.al. [23] which is based on the Lee- Wick model 



of effective QCD. They calculated the Lee- Wick potential at finite temperature to 
obtain the following nucleation rate 

lit) = T 4 J exp (-S/T) (14) 



with 

S 



P = ^-(T 4 -T 4 ) (15) 

where a = lOOMeV and m = 939MeV. 

Csernai and Kapusta [24|] proposed a nucleation rate which is of the form 



2tt ( mal \ 3 1 
7tt 2 ■ 



I(t) = r T exp[-A/(l-t/t c ) 2 \ (16) 

with 

' 16 ™ 3 (17) 



A8B 2 T C 

where a is the surface tension and r T is a temperature dependent constant. They 
have calculated the pre-exponential factor from an effective field theory of QCD. 
This nucleation rate is the same as the general form proposed by Landau and Lif- 



shitz [E5], apart from the pre-exponential factor. 
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a 

(MeV/fm- 2 ) 


T c 
(MeV) 


minimum radius 
(meter) 


N QN 


Cottingham et. al. 


100 


1.66 


7420 




150 


0.083 


1.7 x 10 7 


10 

Csernai et. al. 


100 


0.117 


2.1 x 10 7 




150 


0.0096 


3.8 x 10 10 


50 

Csernai et. al. 


100 


1.25 


1.7 x 10 4 




150 


0.0882 


1.4 x 10 7 



Table 1: Different values of minimum radius and A?£ for the different nucleation 
rates. 
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Let us discuss the results obtained so far. In table I we have shown the de- 
pendence of the minimum radius of a quark nugget and the number of QNs within 
the horizon just after the QCD phase transition on the value of T c for different nu- 
cleation rates. For the nucleation rate proposed by Csernai and Kapusta, we have 
varied the surface tension from a = 10 — 50MeVfm~ 2 . We have found that the 
minimum radius varies from 9.6 x 10~ 3 meter to 1.66 meter. 

In fig. |l] we have plotted the distribution of QN, as a function of ub using 

the nucleation rate proposed by Cottingham et.ai, for different values of T c , where 
Hb is the baryon number content of a single QN. We see that for T c = lOOMeV^ 
distribution of QN peaks at baryon number ~ 7 x 10 45 and there is almost no QN 
with baryon number larger than 10 47 . For T c = 150MeV these values are 10 42 and 
10 43 respectively. In figs. similar results have been shown for nucleation rate 
proposed by Csernai and Kapusta, for a = 10 and 50MeV fm~ 2 , respectively. We 
see that for a fixed T c , the minimum radius increases with increase in a. Figure 

3 shows that the maximum number of nuggets are around a baryon number of 
~ 2.5 x 10 42 at T c = lOOMeV whereas the number of nuggets goes to zero after 

4 x 10 43 . Similar results for T c = 150MeV are also shown in the figures. We should 
mention at this stage that the lower cutoffs that we have obtained here are certainly 
allowed by the study of Madsen et. al. jl2] within the reasonable set of parameters. 

The question which arises next is whether these nuggets will survive till the 



present epoch. Earlier studies [17] have shown that the nuggets having baryon 
number less than 10 42 will not survive till the present epoch. This suggests that 
all the cases considered here for T c = lOOMeV will give stable nuggets. However, 
from figure 1 and 3, it can be seen that, for T c = lbOMeV, some of the nuggets will 
not survive. Fig. 2. suggests that for a = 10MeVfm~ 2 , none of the nuggets will 
survive when T c = lbOMeV. Given the present state of the art, there is no way to 
choose any one of the possibilities as the preferred one. We should therefore con- 
sider the situation that while some nuggets may indeed be stable and constitute cold 
dark matter, some smaller nuggets may evaporate, creating sizeable baryon inhomo- 
geneities. In the next section we will study the evolution of these inhomogeneities 
with time/temperature. 
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3 Evolution of baryon inhomogeneities due to evap- 
orated quark nuggets 



Our aim in this section is to study the implications of those QNs that do evaporate 
away, assuming, of course, that they were formed in the early universe. When a 
QN dissociates into nucleons, the latter initially form a clump with high baryon 
overdensity relative to the density of baryons in the ambient universe. The baryon 
density in the clump then gradually decreases as various physical processes tend 
to 'flatten' the clump. We study the evolution of the highly non-linear baryonic 
inhomogeneities represented by these high density clumps due to dissociated QNs 
created after the epoch of quark-hadron phase transition at T ~ lOOMeV, T being 
the temperature of the universe. 

The evolution of large, non-linear baryon inhomogeneities in the early uni- 
verse has been studied in detail recently, especially in the context of possible cre- 



ation through electroweak baryogenesis process [^6|, |27|, of large baryon inhomo- 
geneities during the epoch of a possible first-order electroweak symmetry-breaking 
phase transition at T ~ lOOCreV. The single most dominant physical process that 
determines the evolution of large baryon inhomogeneities in the early universe before 
the epoch of neutrino decoupling (at T ~ lMeV) is the so-called "neutrino infla- 
tion". Any large baryonic clump in pressure equilibrium with the ambient universe 
would have a slightly lower temperature inside the clump relative to the temper- 
ature of the ambient surroundings, due to the excess pressure contributed by the 
excess baryons inside the clump. As a result, heat would be conducted into the 
clump from the ambient medium. The particles most efficient in conducting heat 
into the clump are the neutrinos which have, by far, the largest mean free path 
(MFP) amongst all the relevant elementary particles. As neutrinos cross the clump 
they deposit energy into the clump thereby heating up the clump. The clump then 
expands in order to achieve pressure equilibrium under this changed circumstance, 
and so the baryon density within the clump decreases as the clump expands. This 
process of expansion ("inflation") of the clump due to neutrino heat conduction 
continues until the neutrinos decouple at around T ~ lMeV. For any given size of 
a clump, the time scale on which a clump achieves pressure equilibrium with the 
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surroundings is essentially the hydrodynamic expansion time scale or the time taken 
by sound to traverse the clump, which can be shown to be smaller than the heat 
transport time scale for neutrinos. It is, therefore, a good approximation to treat the 
evolution of the clump as going through a succession of pressure equilibrium stages 
with decreasing density inside the clump. After neutrino decoupling the evolution 
of the clump is determined mainly by the process of baryons slowly diffusing out of 
the high-density clump to the ambient medium. 

In this section, we study the evolution of the baryonic clumps created by evap- 
orated QNs under neutrino inflation. It is to be mentioned here that the linear re- 
lationship, assumed in ref.[^C], between the baryon overdensity within a clump and 
the fractional temperature difference of the clump relative to the ambient tempera- 
ture, turns out to be invalid in our case of extremely large initial baryon overdensity 
created by the evaporated QNs, as we discuss below. As a consequence, we need to 
numerically solve the full non-linear pressure equilibrium equation for a clump in 
order to obtain the relationship between those quantities. Furthermore, the initial 
baryon overdensity within the clump in our case can be so large (e.g. ~ 10 12 ) pT 



that baryon-to-entropy ratio within the clump could be initially greater than unity 
in which case the dominant contribution to the MFP of neutrinos would come ini- 
tially from neutrino-nucleon scattering rather than from neutrino-lepton scattering 



assumed in ref.|26|. The above two considerations make a straightforward applica- 
tion of the results of ref.|26] invalid in our case of large baryonic inhomogeneities 
due to evaporating QNs; hence the need to do an ah initio calculation for inhomo- 



geneities with initial overdensity significantly larger than those studied in ref.[^| 
In this respect, we believe the calculations in this paper, although done in the spe- 
cific context of inhomogeneities due to quark nuggets, have much wider validity and 
application. We would like to mention at this stage that we are interested to study 
only the neutrino inflation process which will be operative till T = lMeV. So we 
restrict our calculation upto that temperature. Also, the only difference from the 



work of ref. |26j is that we have solved the full non-linear pressure equation. 

The pressure equilibrium equation for a baryonic clump with baryon number 
density p* B = 8nPb and temperature T* = T(l — 5T) in the background universe at 
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temperature T and baryon number density ps can be written as 

Pb T* + p eff (T*)aT* 4 = p B T + ^ e// (T)aT 4 (18) 

where g e ff(T) is the effective number of relativistic degrees of freedom in the universe 
at temperature T contributing to the energy density and pressure, and a = 7r 2 /30. 
The baryons within the clump as well as outside are assumed to be ideal gases 
of non-relativistic particles with pressures p* B T* and p B T, respectively. Assuming 
Sn >1,6T< 1, and g eff (T*) « g eff (T), we get from eq.( [Lj), 

ST - (19) 
1 + V6 N 

where rj = Pb/s, s being the entropy density. The baryon-to-entropy ratio rj in 
the universe is essentially constant for the temperature range of our interest, s ~ 
2.6 x 10 8 Cl B /L h~ 2 , where is the baryonic mass density in the universe in units of 
the closure density, and h = H / (100km, sec' 1 , mpc' 1 ) , H being the present value 
of the Hubble constant. 

Eq.(|P9|) shows that if the baryon overdensity 5n in the clump satisfies the con- 
dition t]6n <C 1, then 5T ~ t]5n, i.e., 5T is linearly proportional to 5n On the 



other hand, for overdensities satisfying 7)5^ ^> 1, Eq.([TP|) gives 6T ~ 1, which is 
inconsistent with the assumption 5T 1 under which Eq.(|19"D is derived. Clearly, 
then, for sufficiently large overdensities for which 5n > T] -1 , the assumption ST 1 
is not valid, and so we need to solve the full non-linear pressure equilibrium equa- 
tion, Eq . (|18D , to obtain the relationship between 5T and 5^. This is the essential 



difference between our work and that of ref.|26|. The result is demonstrated in figure 



Now, for a given overdensity 5^ of the clump at some time t when the temper- 
ature of the universe is T, the rate of energy deposited into the clump by neutrinos 
depends upon whether the size L (= 2R, R being the radius) of the clump (assumed 
spherical) is larger or smaller than the MFP (A^), of neutrinos through the clump 
at that time. For X u < L the clump will not inflate by any significant amount be- 
cause the energy deposition by ambient neutrinos will occur mainly in a thin surface 
layer of the clump leaving the bulk of the clump unaffected. Indeed in this case the 



heating of the clump will be governed by slow diffusion |27| of neutrinos inside the 
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clump. However, X v in the early universe increases rapidly as T decreases, typically, 
X v oc T~ 5 . This means that a clump of any given size L will quickly come within 
the "neutrino horizon", such that X u becomes larger than L before any significant 
neutrino inflation of the clump takes place. Indeed, most of the neutrino inflation 
of the clump will take place when L ~ X u and L > A„. The evolution of 5n with 



time (temperature) is governed by the following differential equations 0, ^£ 



dt R p 

for L ~ X u , 



5T5 N (20) 



= -7^--^ (2i) 

dt 4 \ v p y ' 

for L > X u . 

The typical values of the overdensity 5n and size L of those overdensities ex- 
pected from QN evaporation, are calculated by Iso et al [PI . The values of Sn 
could be as large as 10 12 and R ~ 10 cm. Since nothing is known about the initial 
overdensity 5^ and the length scale L, we study the evolution of the baryon over- 



densities for various initial values of 5n and L by solving equation (gy) and (|2l|). 
In figure |] we have shown the importance of considering the non-linear term. It 
can be seen from the figure that at high 5T/T the linear relation breaks down quite 
substantially. The nuggets which will not be stable against evaporation will form 
highly dense baryonic lumps. We have studied the evolution of these lumps with 
time. Neutrinos play an important role in the evolution of these lumps up to 1 MeV. 
The results are shown in fig. [5[ Lumps with initial overdensity < 10 8 is not affected 
by the neutrino conduction. The final values of the overdensity is smaller for higher 
initial temperature. 

From the above discussion and fig. 5 it is clear that some overdensity is left 
out after the neutrino inflation which is of the order of 10 7 . This overdensity, as 
it looks, is a sizeable amount. The baryon diffusion starts dominating after the 
neutrinos fall out of equilibrium (T ~ lMeV). From T = lMeV to the beginning 
of the nucleosynthesis i.e. T = O.lMeV baryon diffusion is the most dominant 
process as far as the dissipation of the overdensities are concerned. If the baryon 
diffusion lengths are larger than the typical size of the inhomogeinities then the 
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overdensities will be washed out due to this process. As a result these objects will 
not alter the standard big bang nucleosynthesis scenario. This has been shown 
for the baryon inhomogeinity created at the Electro- Weak scale by Brandenberger 
et. al. |30fl . These findings once again supports the existence of nuggets. If the 
evaporating nuggets would have left very high asymmetries in the universe then the 
observed He 4 abundance, which is thought be very well determined, would have 
been violated, a scenario not very comfortable with the survival of quark nuggets. 

4 Conclusion 

In this work we have estimated the abundance of quark nuggets in various nucle- 
ation scenarios with different values of critical temperature and surface tension of 
the bubble. We have found that within a reasonable set of parameters QNs may 
be a possible candidate for cosmological dark matter. The evolution of baryon 
inhomogeneities, formed due to the unstable QNs have also been studied. 



5 Appendix 



In this appendix we follow Ref. [19] to solve the following integral equation, 



poo A-jr 

P(X;t)= ( v - X ) 3 F( m t)dr) (22) 

Ja+X 3 

Differentiating the above equation for four times we get 

^P^(X) = -D(ad x ) F(X) (23) 

47T 

where 

<9 3 f) 2 8 

D(ad x ) = a 3 — - 3a 2 — - 6a— - 6 (24) 

Let 

C(F(X)) = F(p) = / F(X)e- pX dp (25) 
Jo 

where C(A) is the Laplace transform of A. Now using the Laplace transformations 
of the derivatives we have 

Ap(4)( X ) = D(ap)F(p)e- ap (26) 

4:71 
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where we have neglected an arbitrary constant. Now, 



F(X) 



3 l rc+ioo p(*)(p) 



4tt 2m Jc-ioo 
3 r°° 



■ p(x - a) dp 



4:71 J0 ' y " 2Tli Jc 



D(ap) 

I i-c+ioo P p{X-a-T)) 



D(ap) 



dp 



-PV>(X -a)- -P(X -a) 
a 

+ — d V P{ V + X -a){\e-^ /a 
cr Jo ^ 



(27) 



F(X) 



3 6{X - a) 
Ana 3 



-P\X - a) 



3P(X - a) 



a 



/>00 . . 

+ — / dr]P(ri + X-a) f Ae ( - A??/Q) + + uoe^ 1 ^ ) 

cr Jo L > 

where A, u and uj are the solutions of the equation 



x s - 3x 2 + 6x - 6 = 



In terms of the variable z Eq. (^) looks like 



(28) 



(29) 



F{z) 



3 6(z - a)R{ti 
4-KoPvHi 



-P'(X - a) 



3P(X - a) 



a 



1 f°° r 
+ — / dr]P(r] + X - a) Ae 
cr Jo 1 

R(t M 



(-\r)/a) _|_ we (-a)))/a) _|_ ^(-urj/a) 



with 



at>t c 



R(t t ) 



(30) 



(31) 
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Figure 1: Distribution of QN, /(%), as a function of n B using nucleation rate 
proposed by Cottingham et. al. 
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Figure 4: The relation between baryon inhomogeneity and temperature difference. 



The dotted line corresponds to the solution of eq. (18) and the solid line corresponds 
to the linear approximation as discussed in the text. 
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Figure 5: The evolution of baryon inhomogeneity with temperature, for different 
sizes of the clump and different initial temperature, due to neutrino inflation. 



